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Abstract — We consider the pulse design problem in multicar- 
rier transmission where the pulse shapes are adapted to the 
second order statistics of the WSSUS channel. Even though the 
problem has been addressed by many authors analytical insights 
are rather limited. First we show that the problem is equivalent 
to the pure state channel fidelity in quantum information theory. 
Next we present a new approach where the original optimization 
functional is related to an eigenvalue problem for a pseudo 
differential operator by utilizing unitary representations of the 
Weyl-Heisenberg group. A local approximation of the operator 
for underspread channels is derived which implicitly covers the 
concepts of pulse scaling and optimal phase space displacement. 
The problem is reformulated as a differential equation and the 
optimal pulses occur as eigenstates of the harmonic oscillator 
Hamiltonian. Furthermore this operator-algebraic approach is 
extended to provide exact solutions for different classes of 
scattering environments. 

I. Introduction 

Pulse shaping in multicarrier transmission is a key ingredi- 
ent for high rate wireless links. Furthermore it is the standard 
tool to mitigate the interference caused by doubly dispersive 
channels. Most multicarrier schemes like conventional OFDM 
exploiting guard regions (a cyclic prefix), pulse shaped OFDM 
and OFDM/OQAM can be jointly formulated. Hence we focus 
on a transmit baseband signal s(t) given as 

(1) 



s{t) 



• where i is the imaginary unit and 7„ 



def 



\nT,raF)l are 



X 
S3 



time-frequency shifted versions of the transmit pulse 7, i.e. 
shifted according to the lattice TZ x FZ. It is also ben- 
eficial to consider different lattice structures [1] on which 
our contribution will apply as well. The time-frequency (or 
phase space) shift operator Si^r.u) is intimately connected to 
unitary representations of the Weyl-Heisenberg group as we 
will elaborate later on. Therefore (1) is also known as Weyl- 
Heisenberg or Gabor signaling. 

The coefficients Xmn in (1) are the complex data symbols 
at time instant n and subcarrier index m with the prop- 
erty E{a;a;*} = I (•* means conjugate transpose) where 



y . The indices {ran) range over some 



doubly-countable index set I, referring to the data burst to 
be transmitted. We will denote the linear time-variant channel 
by Ti. and the additive white Gaussian noise process (AWGN) 
by n{t). The received signal is then 



with S(t, v) being a realization of the "channel spreading 
function". In practice S(t, v) is causal and has finite support. 
We used here the notion of the WSSUS channel. In the 
WSSUS assumption the channel is characterized by the second 
order statistics of S(t, v), i.e. 



r{t) = {■Hs){t) + n{t) = 



S(t, v){S(r,u)s){i)dTdv + n{t) 

(2) 



E{S(t, i^)S(r', v')} = C(t, v)5{t - t')6{v - v') 

where C (r, v) is the scattering function. Without loss of 
generality we assume ||C||i = 1. To obtain the data symbol 
Xki the receiver does the projection on gki = S(j_T,kF)g^ i-S- 

Xki = {gkur) = j g^i{t)r{t)dt 

def 

By introducing the elements HkUmn = {gku'H'^mn) of the 
channel matrix H G C-^^-^, the multicarrier transmission can 
be formulated as the linear equation x — Hx + n, where n 
is the vector of the projected noise having a power of per 
component. We assume that the receiver has perfect channel 
knowledge (given by S(t, v)), i.e. single carrier based equal- 
ization in the absence of noise would be = Xki/Hki^ki, 
with 

Hum = (gkhTi-lki) = j j T,{T,v){gkuS(r,u)lki)dTdv 

S(t, i^)e^'2'^(^'^-^-'''^) Ag^(T, v)dTdv 

where Ag^{T,v) ~ {g, S(^t-m)1) is the cross ambiguity func- 
tion of the pulse pair {3,7}. 

II. Problem Statement 
Considering only single carrier equalization, it is natural to 
require a '= \HM,ki\'^ ('^he channel gain) to be maximal and the 
interference power h = Yj{ki)^(mn) Wki.mn? to be minimal 
as possible. This addresses the concept of pulse shaping. 
However to be practicable, the pulses should be adapted to 
the second order statistics only, given by C (t, v) and not 
to a particular channel realization H{t,v). Hence, we aim at 
maximization of 

E^{a} 



SINR 



def 



r2 + E„{6} 

by proper design of 7 and g. Up to very few special cases 
the analytical solution of this global optimization problem 
(jointly non-convex in (7, g)) is unknown. However numerical 
optimization methods are presented in [2], [3], [4]. Following 
our previous work [4] we simplify the problem by proposing a 



relaxation, which separates the problem into two steps. Upper 
bounding Ei-j^lb} < — £7^(0} gives a lower bound on 
SINR (see [4]), where is the so called Bessel bound of 
{7mn} [5]. In this paper we focus on the first step only 
where E7^{a} should be maximized. This gives the following 
optimization problem 



arg max 

|7||2 = ||S||2 = 1 



(3) 



def 



where d/i = C(r, v)dTdv. In this context it was first intro- 
duced in [6] respectively [7], but similar problems already 
occurred in radar literature much earlier In particular for the 
elliptical symmetry of C(t, v) Hermite functions establish 
local extremal points as found in [7]. The scaling rule for fixed 
pulses was studied in [7], [8]. Also it is possible to find a close 
relation to the channel fidelity and minimum output entropy 
states in quantum information theory as we will show later 
on. In particular the important class of Gaussian scattering 
profiles (corresponding to classical bosonic quantum channels) 
was already addressed in [9] and [10]. 

Out of the scope of this paper is the second step, in which 
the minimization of B^ (which depends on ^(°p')) is achieved. 
This well known procedure [1] (in the case of TF > 1 that 
is to find the "nearest" orthogonal Gabor basis with respect 
to the £2-norm) is also described in [4]. Unfortunately the 
resulting pulses will be in general again a suboptimal solution 
of (3). See [11] for a discussion of this problem. Nevertheless, 
this separation and therefore (3) opens up analytical insights 
into the pulse design problem. 

III. Contributions 

The original formulation of the pulse design problem in 
(3) hides the internal group structure induced by the time- 
frequency shift operators. In this paper we derive a lower 
bound for the optimization functional (3) on which we can 
exploit this structure explicitely. Moreover we sketch that the 
results will hold in the direct problem with minor restrictions. 
We present an operator-algebraic reformulation by utilizing 
representation theory of the Weyl-Heisenberg group. Our ap- 
proach relates the optimal pulses to approximate eigenstates of 
pseudo differential operators. The procedure naturally embeds 
the concepts of pulse scaling and optimal time-frequency 
offsets (or phase space displacement). Then we extent our 
framework to provide exact solutions for the class of Gaussian 
scattering profiles. Because the underlying theory is partially 
not very common in multicarrier community we will give a 
short introduction to the few properties we will need for our 
investigation. More details can be found in [12]. 

A. The Weyl-Heisenberg Group and Pseudo differential Oper- 
ators 

The two families of shift operators S(r,o) ™d S(^o,u) are 
unitary representations of the group corresponding to the real 
line R with addition as group operation. The extension to 
in the sense of 




is not closed because of the phase factor Closeness is achieved 
by introducing the torus (T) as the third variable, i.e. 

e *(a,,3)-e - e ^ *(a+7,/3+5) 

The corresponding group H = M x M x T with the group 
law [a, P, (f)){"i, (5, V') = (a + 7, /3 + (5, + ?/> — a5) is called 
the (reduced') polarized Heisenberg group (HG). The HG can 
be represented as a group of upper triangular matrices by the 
group homomorphism 



where the group action is matrix multiplication. The matrices 

h{a,(3,4>) = H{a,(3,4>) - 1 written with d = (1,0,0), 
X = (0,1,0) and e = (0,0,1) as h{a,l3,(j)) = ah{d) + 
l3h{x) + (j)h{e) are clearly isomorphic to and with the 
matrix commutator they turn into a Lie algebra. The Lie 

def 

bracket in this case is [(a, /?, 0), (7, J, -0)] = (0,0, — 
(Sj). Due to the bilinearity of the Lie bracket this can 
be shortly written as the Heisenberg Commutation Rela- 
tions, i.e. [d, x] ~ e [x, e] = [d, e] = 0. That this is exactly 
the Heisenberg algebra connected to the HG follows from 
h{a,P,(l)f = /i(0,0,a/3) and h{a,l3,(t>Y = for n > 2. 
The exponential map of the matrix h{a, (3, <f>) is then given as 



= f: M^iM: = 1 + ^) + hiO, 0, a/3) 



(5) 



z27rQ(5 t 



'(a+7,;9+<5) 



(4) 



Thus, it maps the Heisenberg algebra to the unpo- 
larized HG. The series expansion is finite (the ele- 
ments h{a,(3,(j>) are nilpotent endomorphisms). Returning 
to the polarized Heisenberg group we transform finally 

(a, /?, (p) = H{0, 0, -ia/3)e''("''5''^). 

To establish the connection to 8(^^,0) considered as op- 
erators on iS(R) (the Schwartz space of rapidly decreasing 
functions) we have to switch to the so called Schrddinger 
representation. In this picture the hermitian operators X and 
D with 

2m 

setup a basis representation for the Heisenberg Lie al- 
gebra. The skew-hermitian operators 2TiiX (generates the 
frequency shifts), 2-KiD (generates the time shifts) and 
2TTiE = 2TTi (E is the identity) correspond to x,d and 
e. They give again the Heisenberg commutation rules, 
hence linear combinations of them fulfill the same Lie 
bracket (the commutator of linear operators) and consequently 
(r, V, s) dp{T, v, s) = 2TTi{s + vX + tD) is again a Lie 
algebra isomorphism for the Heisenberg algebra. As in (5) the 

'The addition in third component is taken to be mod 1. Otherwise this 
yields the (full) polarized Heisenberg group with non-compact center. 



HG is then given by exponentiation, i.e. the so called Weyl 
operator is given as 



With p{t, v) = p{t, ly, 0) we have S^^.u) = e'^'^^pi^r, v) = 
^ititu^2-ki{vX-tD) ^ j g integrals over shift operators as in (2) 
are in fact pseudo differential operators [12] of the following 
spreading representation (the Weyl transform) 



cr(D,X) 



(7(t, i')e 



2iri(yX+TD] 



drdv 



a(T^ v) is called the spreading function (or representing func- 
tion, i.e. the 2D Fourier transform of the symbol a[d,x) of 
the operator ct(JD, X)). 

B. The WSSVS Pulse Design Problem 

Straight forward calculation shows now that the squared 
magnitude of the cross ambiguity function |Ag^(r, can 
be written in the following form 



(6) 



where G (F) is the (rank-one) orthogonal projector onto g 
(7). By that transformation we emphasize that F undergoes 
a linear transformations before being projected onto g. This 
special kind of linear transformation is also called a unitary 
evolution, which preserve the spectrum of F (in our case the 
rank). This obviously does not hold in generality if taking the 
sum over different unitary evolution of the same argument. 
Hence, we collect them together by defining affine maps A 
and A such that 



def , 



(7) 



The main reason for this reformulation is the notion of 
completely positive maps (CP-maps) [13] which directly 
apply on the pulse design problem. CP-maps like A{-) 
received much attention due to its application in quan- 
tum information theory. Before going more in detail, let 
us define 71 as the set of trace class operators. The set 

def 

Ml = {z \ z & 7i, z = z* , z > 0,Tr z — 1} is a convex sub- 
set of Ti. With Z we will denote the extremal boundary of 
Ml, which is the set of all orthogonal rank-one projectors. 
With the definition of A in (7) follows that A is adjoint of A 
with respect to the inner product TrX*Y. Due to ||C||i = 1 
both maps are trace preserving 'TrA{X) = TrX. Moreover 
they are hermiticity preserving A{X)* = A{X*) and entropy 
increasing X >- A{X) [y is the partial order due to eigenvalue 
majorization). The complete positivity and the trace-preserving 
property is ensured by 



With this framework we can write now the optimization 
problem as 

max TrGAiT) ^o^ 



where F represent the transmitter and the CP-map A{-) 
represent the "averaged" action of the channel and G is the 
receiver. This formulation is similar to the channel fidelity 
in quantum information processing. In fact - the problems 
are equivalent if considering so called pure states. The initial 
preparation of a pure quantum state (the symbol to transmit) 
is represented by a so called rank-one density operator (in our 
case F). The quantum channel is represented by a CP-map A{-) 
having again a density operator as its output. The measurement 
(the detection of the transmitted symbol) is performed in our 
case with G. Obviously either G or F can be dropped in the 
optimization, i.e. 

max TrGyl(F) = max ||A(F)||oo = max j| A(G)||oo (9) 

where ||-|joo denotes the operator norm. This measure rep- 
resents the maximum achievable purity of the output of a 
quantum channel with pure states as input. Turning back to the 
language of WSSUS signaling, this represents the maximum 
achievable "energy" which can be collected by a single pulse if 
communicating with the optimal pulse 7 over a large ensemble 
of WSSUS channels. CP-maps over the Heisenberg group 
have some more important properties. One is the covariance 
property with respect to group elements which follows from 
(4), i.e. 



(10) 



The physical meaning is that (8) is invariant with respect 
to common time-frequency shifts of G and F. A trivial 
but important conclusion is that Weyl-Heisenberg (Gabor) 
signaling is a reasonable scheme, which guarantees the same 
performance on all lattice points. Alternatively it can be 
viewed in the quantum picture as symbol alphabet of pure 
states achieving all the same fidelity. From (10) follows 
furthermore that different maps Ai and A2 commute, i.e. 
Aio A2 ^ A20 Ai. Coming back to the formulation of pulse 
design problem in (9) we can finally relax the constraint set 
from Z to Ml which gives 



max |1A(F)|| 
reMi 



max 

GeMi 



\A{G)\\. 



(11) 



provided by the convexity of ||-||oo and linearity of A{-). To 
the authors knowledge this reformulation of the pulse design 
criterion as a convex maximization problem seems to be new. 
Without further investigations of the analytical structure of 
A{-) such global-type optimization problems are in general 
difficult to solve. Therefore we will emphasize in the following 
more on the Heisenberg group structure contained in A{-). 

C. The Schrodinger Representation 

The connection between Weyl operators (the unitary repre- 
sentations of the Weyl-Heisenberg group in the Schrodinger 
picture) and S{^r,u) will reveal the fundamental role of Gaus- 
sians in WSSUS signaling. We will show this first in a 
simpler lower bound analysis which mainly admits the same 
maximizer as the original problem (given in the appendix). 
Thus, coming back now to (3) and let (tq, j^o) be an arbitrary 



offset between g and 7 in the time-frequency plane, hence we 
define 7 = S^ro.vo)!- 

E-H{a}=y" \Agy{T,iy)\'^dfi^ J \{g,S(r^ro,u-uo)l)\^d^i 



>{j I(.9.p(-t + to,i^-i^o)7>Mm)^ 



= I (.9, [ j Pi~r + T,,v^ v,)dy]^) p I (.9, £7) ? 

In the latter we used Jensen's inequahty^ (J d/i = ||C||i = 1, 
see also [4]). We will use now (12) for further analytical 
studies. The bound becomes sharp iff S^Ag^{T,v) G M is 
constant on supp C for some ^ G T, hence is well suited for 
underspread channels. The operator £ is a pseudo differential 
operator with spreading function ^(r, v) = — C(ro— t, v+vq). 



C 



Local approximation: The nilpotent property with respect to 
the matrix product celebrated in (5) unfortunately does not 
translate into the Schrodinger picture, so that 

%,.) = e'^^"V(-T,i^) 

= e^'^'^ll + dp{-T, v) + idp(-T, vf] + o(2) (13) 

def 

with the hermitian operator K = vX — tD holds only as 
an approximation (for t and v being small), i.e. gives a local 
approximation of C which is 

L = Coo + 27rz(CoiX-Cior>) 

- 27r2(Co2X2 + C20-D' - Cii[XD + DX]) 

where Cmn = // C{t^ v){t — tq)'"^ {v — vq)"" are the moments 
of the scattering function around (roji'o). Because X and D 
are hermitian operators, L is hermitian too if C„„iz'"+" e M 
for m,n = 0,1,2. In this case the optimization problem is 
an eigenvalue problem. Moreover then it follows that g = 
aL^ for some a S C, because only in this case equality in 
|(5i,L7)| < ||g|l2||i7||2 is achieved. L can be made hermitian 
if we choose tq = 1|tC||i and i/q = ||i^Cj|i, so that Cio = 
Coi = 0. Thus we have 

L = Coo - 27r2[Co2X2 + €2^0^ - Cii{XD + DX)] 

which is an hermitian differential operator of second order. 
With 

(d^f)it)'^— fit/a) 

def def 

we define now dilated functions ga = dag, = daj and the 

def 

dilated operator La = daLdi/a- Using furthermore that 



daXdi/a — — X 

' a 

daDdi/a ~ ceD 
^It can be shown that Ag-y S Li(/x). 



(14) 



2n^\^X' + C2oa'D' 



we get 

{ga, La^a) = (Sct, C'oO 

-Cii{XD + DX)]^a) 
and with = C02/C20 the phase space symmetric version 

La = 27r^ VCo2C2o{k 



[X' + - CiiiXD + DX)]} 



(12) where the constant is k = 



Cqq 



(with our WSSUS 



assumptions follows also Coo ~ !)• For simplicity let us 
assume that the shifted scattering function is separable yielding 
Cii = 0. In the general case the Cn-term can be removed 
using a proper symplectic transformation (see for example 
[14]). The eigenfunction of the so called sub-Laplacian (or the 
harmonic oscillator Hamiltonian) X^ + D^ are the Hermite 
functions /i„, with {X^ + D'^)hn = ^^h^. Therefore it 
follows that 



Lahn = (Coo - 7rVCo2C2o(2n + l))h„ 

Hence in local approximation the maximization problem is 
solved by ho, i.e. g dj/^/io and 7 = S'f^^l^^^^di/aho 
which are both scaled and proper separated Gaussians (the 
ground state of the harmonic oscillator). This is an impor- 
tant (and expected) result for the pulse design problem in 
WSSUS channels. It includes the concepts of pulse scaling 
(by di/a) and proper phase space displacement (by S'^^ 
as natural operations. However this approximations is only 
valid for C02C20 <C 1 (underspread channel), such that (Coo — 
7rVCo2C2o(2n + 1)) > 0. We obtain the same solutions in 
the original problem if we apply this approximation (see the 
appendix). Next we will derive cases where this approximation 
turns out to be is exact. 

Gaussian scattering functions: Let us assume that after 
performing proper pulse scaling and separation the scattering 
function is given as the symmetric Gaussian C (t, f) = 
ag-f a(r^+i.^) ^jjgj.g Q < a £ R. If Of > 1 the channel 
is underspread. It can be shown that then C essentially self- 
adjoint, hence the maximum in (12) is again achieved by 
eigenfunctions of C. Operators having such spreading func- 
tions are contained in the so called oscillator semigroup [15] 
and for a > I they have the representation [12] 

£ _ g-27i-arcothQ(X^+D^) 

Thus we have that C-K^ g-(2n+i)(arcotha)^^^ j^gj^^g 
the optimum of (12). The special case a = 1 can be included 
by observing that then C{t,v) ^ Ahoha{T, i')- Such pseudo 
differential operators perform simple projections, in this case 
onto the span of ho- Note that for (7(t, i/) = {(p , pir , v)i})) 
follows 

(g,a(r>,X)7) = (a,(.g,p(.,.)7)) 

= ((0,p(.,.)V'),(.9,p(-,-)7)) 



(.9, 0) (■0, 7> 



Thus a{D, X) is a rank one projector (orthogonal in the case 
iP = (f>)if (0,^) = 1. 

Finally we conclude that for underspread channels the Gaus- 
sian pulse shape is an approximate solution of (12) which 



becomes more optimal as the support of C decreases. Further- 
more the solution is exact for a Gaussian scattering function. 
In the quantum channel context the same arguments hold for 
coherent states (phase space translated Gaussians). 

Appendix 

We will sketch now that the results obtained from the lower 
bound analysis will hold with minor restrictions in the direct 
problem. Writing the CP-map A{-) using shift operators gives 

where we introduced again an arbitrary offset [tq.vq) as 

~ def 

already done in (12) for the lower bound analysis, i.e. F = 
^(ra,u,,)^S\roMay Using again (13) gives 

-4(r) ~ y (1 + 27rjA' - 27r2i^2)f (1 - 2mk - 2Tr^K'^)dn 

~ def 

where the self-adjoint operator is now K — {{v — i'o)X — 
(r — To)-D). Hence we get the following approximation on the 
optimization functional 

Tr GA(F) w Tr J (^(1 - 4:n^K^){fG + GT) + iir^KYKG 

+ 2TTiK[r,G]jdfi + oi2) 

If we restrict furthermore g and 7 to be real it can be shown 
that this will become 

TrGA(F)wTr J (^{1 - 4,T:^k^)rG + in'^kTKG 

+ 2^:ik[f,G]Jd^l 

Thus from the latter we can separate the following approxi- 
mated version Ai{-) of the CP-map A{-) 

Ai{T) = J (^(1 - 4:TT^k^)r + in^krk + 2TTi[k, r]^dn 

In the next steps we will perform the integration of the three 
integrands. Using 

J k^Tdn =C2oD^r + Go2X^r - Cn{DXr + XDT) 

J krkdfi ^C2oDrD + Co2xrx - GniDfx + xtd) 

J [k,f]dfi ^Cio[D,t]^Coi[X,t] 
we get finally 

^i(r) = Coo - 4^' {C20D[D, f] + Go2X[X, f]) 

+ 47T^(CnD[X,r]+GnX[D,r])') 

+ 2^^(Clo[^>,f]-Col[X,f]) 

We choose again (to,i'o) such that Cio = Coi = 0. Fur- 
thermore let us assume again for simplicity that the scattering 
function is separable around (tqjI'o) yielding Cn = 0. Then 
we will get 

^i(r) = Coo - 4^' (C2oD[D, f] + Go2X[X, f]) 



Next we apply the same dilation procedure used already for 

~ def ~ 

the lower bound analysis. Let F^ = da^di/^- Using again 
(14) gives 

Ai(F) = Coo - 47rVCo2C2o [d[D, f + f ,]) 

This will result in the same pulse scaling rule as from the lower 
bound analysis. Hence it remaines to show that Gaussians 
are the right pulses to perform the scaling, thus we aim at 
maximization of j|Ai(F)||oo which is 

max|jK/2-Z}[r>,fj-X[X,f JIU = 

max||K/2 - {D^ + X2)f„ + £)f„I5 + Xf„X||oo 

where k = ^ ., ^ . It can be shown that for time- 
frequency symmetric Fq and C follows Tr GXTaX = 
Tr GDTaD=Q. With this restriction remains 

max||K/2 - (D^ + X^)fa\\oo 

which is maximized by F^ being a projection onto the 
eigenspace of + X"^ corresponding to the minimal eigen- 
value which is again ho- For Gaussian scattering functions in 
turn the calculation in [10] suggest that this will hold also for 
A{-). 

References 

[I] Thomas Strohmer and Scott Beaver, "Optimal OFDM Design for Time- 
Frequency Dispersive Channels," IEEE Trans, on Communications, 
vol. 51, pp. 1111-1122, Jul 2003. 

[2] D.Schafhuber, G.Matz and F.Hlawatsch, "Pulse-Shaping OFDM/BFDM 
Systems for Time- Varying Channels: ISI/ICI Analysis, Optimal Pulse 
Design, and Efficient Implementation," Proc. IEEE PIMRC-2002, Lis- 
bon. Portugal, pp. 1012-1016, 2002. 

[3] Phil Schniter, "A New Approach to Multicarrier Pulse Design for 
Doubly-Dispersive Channels," Proc. AUerton Conf. on Communication, 
Control, and Computing, Monticello, Oct 2003. 

[4] Peter Jung and Gerhard Wunder, "Iterative Pulse Shaping for Gabor Sig- 
naling in WSSUS channels," Fiftli IEEE Workshop on Signal Processing 
Advances in Wireless Communications, Lisboa, Portugal, 2004. 

[5] Ole Christensen, An Introduction to Frames and Riesz bases. Birkhauser, 
2003. 

[6] Werner Kozek, "Nonorthogonal Pulseshapes for Multicanier Commu- 
nications in Doubly Dispersive Channels," IEEE Journal on Selected 

Areas in Communications, vol. 16, pp. 1579-1589, Oct 1998. 
[7] Werner Kozek, "Matched Weyl-Heisenberg expansions of nonstationary 

environments," PhD thesis, Vienna University of Technology, 1996. 
[8] K. Liu, T. Kadous and A. M. Sayeed, "Orthogonal Time-Frequency 

Signaling over Doubly Dispersive Channels," IEEE Trans. Info. Theory, 

vol. 50, pp. 2583-2603, Nov 2004. 
[9] V. Giovannetti, S. Guha, S. Lloyd, L. Maccone and J.H. Shapiro, 

"Minimum output entropy of bosonic channels: a conjecture," Phys.Rev, 

vol. A 70, no. 032315, 2004. 
[10] CM. Caves and K. Wodkiewics, "Fidelity of Gaussian Channels," 

arxiv:quant-ph/0409063, Sep 2004. 
[11] T. Strohmer, "Approximation of dual Gabor frames, window decay, 

and wireless communications," Appl.Comp.Hann.Anal., vol. 11, no. 2, 

pp. 243-262, 2001. 
[12] Gerald B. Folland, Harmonic Analysis in Phase Space. Princeton 

University Press, 1989. 
[13] W.F. Stinespring, "Positive functions on C*-algebras," 

Proc.Amer.Math.Soc, vol. 6, pp. 211-216, 1955. 
[14] Micheal E. Taylor, Noncommutative Harmonic Analysis. American 

Mathematical Society, 1986. 
[15] R. Howe, "The Oscillator Semigroup," Proc. Symp. Pure Math., Amer. 

Math. Soc, vol. 48, pp. 61-132, 1988. 



